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ABSTRACT 

The  volume  of  the  research  effort  currently  under  way  to  characterize 
the  constitution  of  thermoelastic-plast ic-creep  materials  is  ample  evidence 
of  its  importance.  Although  significant  work  has  been  done  to  obtain  working 
constitutive  models,  in  many  cases  the  theory  has  not  been  cast  with  the  re¬ 
maining  field  equations  into  a  working  continuum  model.  Tn  this  paper  the 
author  will  consider  models  from  classical  plasticity  and  nonlinear  visco¬ 
elasticity  theory  for  applicability  to  the  conservation  of  momentum.  Based 
on  physical  accuracy  as  well  as  computational  efficiency  the  author  will  assess 
the  feasibility  of  further  development  of  each  model. 
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INTRODUCTION 

Current  growth  in  computer  capability  has  brought  a  parallel  expansion 
in  the  research  area  dealing  with  the  predicted  response  of  materials  whose 
constitution  is  nonlinear.  One  result  of  this  tremendous  explosion  of  tech¬ 
nology  has  been  the  need  for  increasing  specialization  among  scientists  and 
engineers.  This  specialization  has  caused  the  division  of  constitutive 
theorists  into  two  distinct  and  often  noninteracting  groups:  the  materials 
scientists  and  the  mechanics  engineers.  Each  group  employs  a  different 
approach.  Whereas  the  materials  scientist  develops  a  theory  based  on  micro¬ 
structure,  the  engineering  mechanist  tends  to  formulate  a  model  from  continuum 
theory.  The  result  is  that  many  micro structurally  based  theories  are  compu¬ 
tationally  weak,  while  continuum  theories  fail  to  satisfy  microstructural 
constraints. 

A  sound  constitutive  theory  must  be  both  physically  and  computationally 
effective.  In  order  to  establish  both  of  these  criteria  it  is  necessary  to 
understand  the  problem  at  hand.  Therefore,  in  the  first  section  of  this 
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paper  the  author  will  define  the  constraints  precisely.  Following  this  two 
models  for  analyzing  thermoelastic-plastic  solids  will  be  proposed.  Finally, 
the  author  will  consider  strengths  and  weaknesses  of  each  in  order  to  assess 
which  model  is  more  suitable  for  further  research. 


PHYSICAL  AND  COMPUTATIONAL  REQUIREMENTS  OF  THE  MODEL 


There  are  four  requirements  which  a  useful  constitutive  theory  must 
satisfy:  thermodynamics,  physics,  mathematics,  and  experimental  data  require¬ 

ments.  In  this  section  each  of  these  requirements  will  be  outlined  for  the 
material  models  to  be  considered  herein. 


The  rmodynamic  s 

Thermodynamics  should  be  used  as  the  starting  point  for  any  constitutive 
theory.  From  this  constraint  will  come  a  general  mathematical  framework  for 
the  model. 

Recall  that  the  result  of  the  first  and  second  laws  applied  to  a  general 
system  in  the  neighborhood  of  the  equilibrium  state  is  the  following  N  non¬ 
linear  equations  of  motion  fl,  2,  3]: 


JH  ,  -  o 

3qi  +  bij  dt  “  Qi  ’ 


where  H,  the  Helmholtz  free  energy,  is  given  by: 


H  =  H(qt ,  T)  -  U  -  TS  , 


O) 


(2) 


U  is  defined  by  the  first  law  of  thermodynamics  to  be  internal  energy,  and 

S  and  T  are  defined  by  the  second  law  of  thermodynamics  to  be  entropy  and 

temperature,  respectively.  In  addition,  and  q.  are  generalized  forces  and 

displacements,  respectively.  The  nonlinearity  in  equation  (1)  arises  from 

the  functionality  of  F  in  equation  (2)  and  the  form  of  b  =  b  .  ,  (q,  ,  T)  . 

*  i  1 J  k 

Under  suitably  small  displacements,  Q  may  be  regarded  as  the  infinitesimal 
stress  tensor,  and  qn  as  the  corresponding  strain  tensor,  where  m,n  =  1,2..., 6. 


The  remaining  coordinates  are  called  hidden  coordinates  due  to  the  fact 
that  the  generalized  forces  for  these  indices  are  identically  zero. 

The  specific  mathematical  form  of  the  constitutive  law  will  come 
from  the  form  of  the  Helmholtz  free  energy  suitable  for  modelling  the 
material  under  consideration.  For  the  thermoelastic-plastic  material 
considered  herein,  it  is  assumed  that  the  free  energy  may  be  expanded 
about  any  equilibrium  configuration  in  a  second  order  Taylor  series 
expansion  with  respect  to  hidden  coordinates  similar  to  that  previously 
proposed  [l.  A]: 


H 


H  +dq  +%a'qq 
e  r  r  rs  'r  s 


(3) 


where 


(q^  T)  =  Helmholtz  free  energy  in  state  e, 


dr  =  dr  (V  T)  E  OH/Dqr)e,  and 


a' 

rs 
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rs 


<v T)  ■ *; 


sr 


(?2H/3q  3q  )  , 
r  s  e 


the  subscript  e  denotes  any  equilibrium  state  for  which  .all  =  0,  and 
subscripts  r,  s,  and  t  range  only  over  indices  corresponding  to  hidden 
coordinates.  Substitution  of  equation  (3)  into  equations  (1)  will  yield 
the  following  equations  of  motion: 
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(4) 
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3\ 


3d  dq . 

t  — ■  q  +  a’  q  +  b'  .  =  Q_  , 

Sq^  r  rsns  mj  dt  m 


for  indices  corresponding  to  observed  coordinates,  and 


d 

r 


+  a'  q  +  b ' . 
rs  s  rj 


=  0, 


(5) 


for  indices  corresponding  to  hidden  coordinates,  where  terms  of  second 
order  in  q^  have  been  neglected. 

Specific  constitutive  equations  will  come  from  physical  constraints 
applied  to  equations  (4)  and  (5) . 

Physical  Constraints 

The  thermodynamic  framework  postulated  in  equation  (1)  is  insufficient 
to  describe  a  usable  constitutive  theory.  One  must  in  addition  define  the 
type  of  material  behavior  that  the  theory  should  predict.  The  material  to 
be  considered  herein  is  a  so-called  thermoelastic-plastic-creep  material 
such  as  aluminum  or  steel.  Mathematically,  this  type  of  material  may  be 
represented  by: 


su  '  su  {IW  T)- 


(6) 


where  S..  is  the  stress  tensor,  E  is  the  strain  tensor,  and  {  }  signify 
ij  mn 

history  dependence.  The  functional  dependence  of  the  stress  tensor  on 

the  strain  history  alleviates  the  necessity  to  include  strain  rate  explicitly 

on  the  right  hand  side  of  equations  (6) .  At  this  point  no  restriction  is 


aa-*uS'-;^r.  . '..--.I. 


necessary  regarding  the  mathematical  definitions  applied  to  stress  and  strain. 

Uniaxial  tests  on  thermoelastic-plastic-creep  materials  will  yield  the 
general  macrophenomenological  behavior  shown  in  Figures  (1)  and  (2).  It  is 
seen  from  the  data  that  several  phenomena  which  should  be  modelled  by  the 
constitutive  theory  are:  (1)  the  Bauschinger  effect  for  reverse  loading, 

(2)  the  nonunique  mapping  from  strain  to  stress,  (3)  the  pronounced  effects 
of  temperature  and  strain-rate,  and  (4)  the  general  nonlinear  behavior  of 
the  material.  It  may  also  be  noted  from  the  figures  that  the  entire  load 
history  must  be  known  in  order  to  evaluate  the  state  of  strain,  and  that  the 
effects  of  temperature  are  quite  similar  to  those  caused  by  strain  rate. 

All  of  the  above  phenomena  can  be  modelled  using  equations  (6) . 

The  microphenomenological  behavior  of  the  material  is  quite  complex. 

The  material  behavior  is  characterized  by  an  elastically  reversible  deformation 
component  which  is  due  to  stretching  of  crystalline  bonds,  as  well  as  a  par¬ 
tially  reversible  deformation  due  to  sliding  of  grain  boundaries  and  certain 
complex  phenomena  such  as  dislocation  glide  and  diffusion  of  molecules  within 
grains.  Elevated  temperature  causes  molecular  excitation  with  resulting  bulk 
expansion.  The  Bauschinger  effect  is  primarily  due  to  microscopic  stresses 
called  back  stresses  which  are  built  up  during  deformation. 

The  most  important  physical  restriction  imposed  by  micromechanics  is 
that  there  is  no  physically  observable  mechanical  difference  between  so- 
called  rate  independent  plastic  strain  and  rate  dependent  creep  strain. 

These  phenomena  have  often  been  described  separately  by  macromechanians 
even  though  they  may  be  attributed  to  the  same  processes. 

Incorporating  the  above  physical  phenomena  into  equations  (6)  and  subjecting 
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them  to  the  thermodvnamic  constraints  imposed  by  equations  (4)  and  (5)  will 
result  in  a  phvsicallv  based  constitutive  theory  for  thermoelastic-plastic 
materials . 

Computational  Considerations 

It  is  important  that  the  constitutive  theorv  founded  on  phvsical  grounds 
also  be  computationally  usable  in  solving  continuum  problems  wherein  the  state 
of  stress  and  strain  varv  spaciallv.  This  will  require  that  the  application 
to  kinematics  and  conservation  of  momentum  can  be  performed  efficiently.  To 
perform  this  portion  of  the  model  development  it  is  necessary  to  he  familiar 
with  the  various  methods  for  satisfying  conservation  of  momentum.  Since  a 
detailed  explanation  of  those  methods  currently  in  use  is  beyond  the  scone  of 
this  paper,  only  one  method  will  be  considered  herein. 

The  procedure  chosen  here  is  the  finite  element  method  applied  to  the 
virtual  work  variational  principle.  The  impetus  for  choosing  this  method  is 
that  it  is  already  widely  in  use  in  computer  codes  such  as  AOCIF.  T  [5], 

NONSAP  [6]  and  ADTNA  [7].  Thus,  development  of  a  new  constitutive  theorv 
will  require  onlv  modification  of  a  relatively  small  portion  of  the  code. 

In  order  to  determine  the  constraints  imposed  bv  this  computational 
scheme,  it  is  appropriate  to  review  briefly  the  method  used  in  the  computer 
codes  mentioned  above.  Recall  that  momentum  is  conserved  by  the  theorem  of 
virtual  work  [61,  given  by 


/ 


S  t+AtSF  ,t+AtdV 
i.i  il 


,t+At 


(7) 


where  .S  represents  a  variation,  t+At  denotes  the  time  at  the  end  of  a  load 
increment,  and  W  represents  the  external  work  due  to  surface  tractions  and  body 


9 


forces,  both  inertial  and  gravitational.  The  constitutive  law  and  kinematic 
eauations  are  then  inserted  in  equation  (7)  to  obtain  a  working  variational 
orincinle.  In  order  to  obtain  a  usable  principle,  however,  it  is  necessary 
to  incrementalize  the  constitutive  and  kinematic  eauations  using 


t+At 

H 


At 


(8) 


and 


E..t+At  =  F. .  . t  +  -5U.  At 
i1  i  i  At 


(9) 


Thus,  the  constitutive  model  must  be  input  in  incremental  or  rate  form.  This 
restriction  results  due  to  two  computational  difficulties,  ^irst,  the  solution 
must  be  obtained  in  small  incremental  load  steps  in  order  to  guarantee  conver¬ 
gence  of  the  nonlinear  equations  resulting  from  constitution  and  kinematics. 
Second,  the  fact  that  stress  and  strain  are  not  uniquely  related  in  the  con¬ 
stitutive  theory  requires  that  the  constitutive  theorv  be  written  in  rate  form 
during  single  signed  load  increments  in  order  to  guarantee  a  unioue  solution.* 
Thus,  equations  (6)  are  replaced  bv 


S.  . 

il 


i.l 


(E 


mn 


T,  k. 


k  ) 
n 


(10) 


*Theoret leal ly  this  restriction  mav  be  circumvented  with  proper  state  variables, 
that  is,  equation  (6)  mav  be  replaced  with  (E^n,  E„in ,  T,  kj,...kn), 

where  again  k  are  state  variables,  however,  to  this  author's  knowledge,  no 
theories  of  this  tvpe  are  currently  available  with  sufficient  accuracv  to  model 
the  thermoelastic-plastic  material. 


in 


where  k_  are  state  variables  obtained  from  the  strain  historv. 

Tt  is  to  be  noted  that  the  above  incrementation  allows  for  stress  to 
be  written  as  a  function  of  strain  and  not  vice  versa  due  to  the  form  of  the 
variational  principle  given  in  equation  (7).  This  requires  that  displacement 
conditions  be  specified  on  the  boundarv  of  the  continuum,  and  that  the  compu¬ 
tational  procedure  will  result  in  a  solution  scheme  for  obtaining  unknown 
displacements  as  opposed  to  stresses. 

The  resulting  variational  principle  is  linearized  and  then  discretized 
using,  the  finite  element  method  [6,  7 ]  to  obtain  a  set  of  coupled  linear  dif¬ 
ferential  equations  in  the  displacements  at  nodal  points.  The  nonlinearities 
in  the  variational  principle  are  accounted  for  by  using  a  modified  Newton-Raphson 
iterative  scheme  [7]. 

Tn  summarv,  then,  it  is  necessarv  to  obtain  a  constitutive  law  which 
satisfies  both  the  computational  constraint  given  by  equations  (10)  and  the 
thermodvnamic  restrictions  in  equations  (4)  and  (5),  as  well  as  the  phvsical 
considerations  previously  noted. 

Experimental  Data  Requi remen ts 

The  amount  of  experimental  data  required  to  determine  a  working  constitu¬ 
tive  model  will  significantly  affect  the  usefulness  of  the  model.  Therefore, 
consideration  should  be  given  to  this  final  constraint  in  determining  the  most 
useful  constitutive  theorv. 
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INVESTIGATION  OF  AVAILABLE  THEORIES 


The  available  theories  for  modelling  response  of  thermoplastic  materials 
are  of  three  distinct  types:  (1)  the  microphenomenologlcal  or  equation  of 
state  theories  [8,  9,  10],  (2)  the  classical  plasticity  theories  | 7,  12,  13], 
and  (3)  the  nonlinear  viscoelasticity  theories  [2,  1 3 ,  14,  15].  The  first 
group  is  derived  primarily  from  studies  of  the  microphysical  characteristics 
of  the  material.  Many  so-called  unified  theories,  wherein  creep  and  plasticity 
are  considered  to  be  a  single  inelastic  term,  fall  within  this  group.  The 
plasticity  theories  are  primarily  macrophenomenological  in  nature  and  are 
generally  based  on  the  notion  of  a  yield  surface  in  stress  space.  The  visco¬ 
elastic  theory  derives  mainly  from  thermodynamics,  with  the  constitutive 
relation  usually  written  in  convoluted  form. 

Although  the  equation  of  state  theories  are  receiving  considerable 
attention  at  this  time,  their  extreme  variation  in  formulation  makes  it  an 
impossible  task  to  consider  them  here.  It  is  hoped,  however,  that  valuable 
information  about  the  remaining  two  groups  may  be  obtained  from  micropheiu '•  :- 
nology.  Thus,  the  theories  considered  here  will  be  plasticity  and  nonlinear 
viscoelasticity.  For  convenience  only  uniaxial  theory  is  presented,  although 
three-dimensional  theory  is  analogous. 

Classical  Plasticity  Theory 

The  author  has  previously  proposed  an  extension  to  isothermal  plasticity 
theory  to  account  for  thermomechanical  and  creep  behavior  [11,  lb,  17]  and 
similar  to  theories  proposed  by  Snyder  and  Bathe  [7]  as  well  ns  Yamada  and 
Sakurai  [12],  That  theory  is  first  reviewed  briefly. 

Recall  that  in  plasticity  theory  one  postulates  the  existence  of  a  yield 
surface  In  stress  space,  which  may  be  extended  to  include  temperature  by 
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F  (o  -ci)  =  k7(/dcP,  T) 


(ID 


where  o  is  the  state  of  stress,  a  is  the  center  of  the  yield  surface  in 

stress  and  temperature  space,  and  k  represents  the  yield  surface  size  as  a 

-P 

function  of  the  history  of  uniaxial  plastic  strain,  /dr  ,  and  temperature,  T. 

Inside  the  surface  F  all  action  is  elastic,  whereas  on  this  surface  the 

-P 

element  undergoes  some  permanent  deformation.  The  parameters  a  and  / dc  may 
be  interpreted  as  state  variables  representing  the  dislocation  arrangement 
and  dislocation  density,  respectively.  Both  may  be  conveniently  obtained 
from  an  appropriate  hardening  rule  and  phenomenological  uniaxial  test  data. 
The  yield  function  [ equat ion  ( 1 1 )]  is  supplemented  with  a  flow  rule,  usually 
given  by  the  normality  condition: 


-P 

d, 


(12) 


where  dl  is  a  scalar  to  be  determined  and  the  superscript  P  denotes  a  rate 
independent  instantaneous  component  of  inelastic  strain.  Tn  addition,  it  is 
assumed  that  the  state  of  stress  is  dependent  on  the  elastic  strain  only; 

P  C  T 

o  =F(r-r  -  r  -  r  )  ,  (11) 


where  the  term  in  parentheses  is  the  elastic  recoverable  strain  caused  bv 
bond  stretching,  and  F  is  the  elastic  modulus.  Tn  addition,  the  superscripts 
C  and  T  denote  rate  dependent  creep  and  temperature  dependent  thermal  strains, 
respectively.  Note  that  equation  (11)  satisfies  physical  constraints  only  if 
the  sum  of  plastic  and  creep  strain  gives  the  total  inelastic  and  nonrecoverahl e 
strain. 
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A  combination  of  equations  (11)  through  (13)  in  rate  form  will  result 
in  the  following  incremental  constitutive  law  during  plastic  loading  [16]: 
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3(7  3a  Jo  3a 


(14) 


where  fi '  is  the  slope  of  the  uniaxial  stress-plastic  strain  diagram,  and 
superscripts  t  and  t+At  denote  parameters  at  the  start  of  the  time  step  and 
the  end  of  the  time  step,  respectively.  The  time  dependence  in  the  elastic 
modulus  is  caused  by  a  temperature  change  during  the  load  step.  Note  that  th 
form  of  equation  (14)  imolies  that  all  parameters  on  the  right  hand  side  may 
be  determined  a  priori,  that  is,  without  recourse  to  the  current  stress  incre 
ment.  Note  also  that  the  above  constitutive  law  in  rate  form  is  not  rotation 
ally  invariant  [18,  19].  Therefore,  its  use  is  restricted  to  inf i7iitesimal 
strains  only. 

The  above  constitutive  law  has  been  implemented  to  the  finite  element 
code  ACUTE  I  [5],  along  with  a  nonisothermal  combined  isotropic-kinematic 
hardening  rule  in  order  to  account  for  the  Rauschinger  effect  during  noniso¬ 
thermal  cvclic  loading.  The  model  has  been  shown  to  give  results  comparable 
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to  experiment  for  some  thermomechanical  load  histories  [16,  17].  Two  signifi¬ 
cant  shortcomings  of  the  theory  should  be  pointed  out  herein.  First,  the 
above  theory  was  obtained  without  recourse  to  thermodynamics.  Second,  the 
assumption  that  creep  and  plasticity  may  be  uncoupled  violates  the  observation 
that  they  are  physically  indistinguishable  phenomena. 

Nonlinear  Viscoelasticity  Theory 

The  nonlinear  viscoelasticity  thoery  is  established  directly  from  the 
thermodynamic  equations  of  motion  [equations  (4)  and  (5)].  These  equations 
are  first  specialized  by  assuming  that: 
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3d 
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where  0  =  T  -  T  and  T  is  the  reference  temperature.  It  follows  that 
K  K 

equations  (4)  and  (5)  may  be  rewritten  as: 
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(20) 


Vi;as  =  as(V  T>  >  0  • 


where  arg»  b^..,  anc*  are  constants.  The  resulting  equations  will  be: 


dqi  - 

aijq-1  +  bij  W  =  Qi  +  Bi°  ’ 


(21) 


where 


<Qi  +  aFainan 


0 


(22) 


(23) 


and  the  reduced  time  parameter. 


(24) 


These  equations  are  identical  to  those  obtained  by  Schaperv  with  the  exception 
of  the  thermal  dependence  of  ap  and  a,,.  Schapery  has  shown  that  using  the 
second  order  free  energy  expansion  in  the  equations  of  motion  (21)  is  analogous 
to  a  generalized  Maxwell  model  with  isothermal  elastic  and  thermal  expansion 
coefficients.  This  analogy,  shown  in  Figure  3,  is  also  true  in  eouations  (21) 
except  that  these  coefficients  are  now  temperature  dependent  through  the  terms 
aF  and  as- 

Now  consider  a  single  Maxwell  element,  as  shown  in  Figure  4.  Tf  one 
considers  small  displacements  only,  then  the  generalized  force  Q  mav  he  re¬ 
placed  by  the  uniaxial  stress  n,  and  the  generalized  displacement  q  is  replaced 
with  the  strain  e.  The  differential  equation  relating  stress  to  strain  can  he 
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aF  =  Optq^.T)  ,  a0  =  a0(q„,T)  ,  as  =  ^(q^T) 

Er,  nr,  oiP  =  CONSTANT 

E'r  =  Spring  coefficient  for  the  r,h  element 

a'r  =  Coefficient  of  thermal  conductivity 
for  rth  element 

Yr  =  Dashpot  viscosity  for  rth  element 


Figure  3.  Generalized  Maxwell  Model  with  Temperature  Dependent  Elements 
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ELEMENT  T 


Em  =  Em(T)  =  aFE 

Hm  =nM^J) 3  nao 
o<,E,T\  =  constant 


as  =  as(6,T) 
aF=  aF  (e,T) 
aD=  a0(e,T) 


Figure  4.  Single  Nonlinear  Maxwell  Element 


found  to  be 


,  ,  T  ,  ,  .  dEu  t  _ 

d  e  _  d  e  +  1  da  1 _ M  _a _  a 

dt  dt  _  t+dt  dt  „  t+dt  dt  _  t  na^ 

em  em  em  d 


where  for  convenience  is  used  rather  than  the  equivalent  form  Ea^,  and  c 
denotes  the  strain  in  the  spring  due  to  a  temperature  change.  Now  assume 
that  the  entropy  production  coefficient  may  be  uncoupled  as  follows: 


.-P 

f (a, c ,T)  +  g(a,e,T) 


where  dc  /dt  is  the  uniaxial  plastic  strain  rate  and  the  nonnegativity  property 
imposed  by  thermodynamics  still  holds.  Since  the  uniaxial  plastic  strain  rate 
is  both  stress  and  temperature  dependent,  equation  (26)  may  be  rewritten  in 


the  following  form: 


=  f  +  g 


[(£),$?♦(£),  Sf]  ' 


A  further  modification  gives: 


f  +  g 


/3o\  dT  /acP\  do] 

»tWt  dt  Wt  dtJ 


Substituting  (28)  into  (25)  gives 


de  _  de 
dt  ~  dt 


n  f +  * 


_1 _  da 

t+dt  dt 


1  dEM  at 
dt  ^ 


T  (^?\  (*s\  is!  I 

L  \30  It  Wt  dt  Wt  dtJI 
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Now  define 


(30) 


(31) 


so  that  the  strain  components  in  the  Maxwell  model  are: 


T 

c 


VT  -  V  ■ 


E 

e 


p 

c 


'  A  8  (^r)t  (if)  t  dT  +/*  (^r)t  01,0  • ,nd 

TR  0 


(32) 

(33) 

(34) 


(35) 


where  TR  is  the  reference  temperature  at  which  the  slider  block  undergoes  zero 

strain.  It  can  be  seen  that  the  viscous  terms  have  been  uncoupled  into  a  rate 

P  C 

independent  plastic  term  (e  )  and  a  rate  dependent  creep  term  (e  ') .  Concep¬ 
tually,  this  can  be  described  as  shown  in  Figure  5.  The  plastic  slider  block 
is  seen  to  slide  according  to  the  stress  and  temperature  inputs,  but  indepen¬ 
dently  of  time.  Physically,  this  analog  describes  the  total  dislocation 
movement  in  a  polvcrvstal  as  being  composed  of  two  segments:  a  portion  of 
dislocation  movement  which  occurs  instantaneouslv,  or  at  least  as  fast  as  the 
loading  is  applied,  and  a  portion  such  as  diffusion  which  occurs  in  a  delayed 
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Figure  5 


£t  =  0<m(T"Tr) 

eE  =  °/EM 


ep  =  ep(cr,e1T) 
ec  =  ec(o\e,T,t) 


Specialized  Nonlinear  Maxwell  Element  with  Rate 
Independent  Slider  Block 


2 


fashion.  However  one  defines  these  terms,  the  end  result  must  be  that  the 
sum  of  the  strains  in  the  plastic  slider  block  and  dashpot  must  equal  the  total 
inelastic  strain. 

If  one  considers  strain  as  the  input  or  independent  variable  it  can  be 
seen  that  equation  (29)  is  a  nonlinear  first  order  differential  equation  in 
the  dependent  variable  a.  There  are  several  accented  procedures  for  solving 
the  equation  for  the  stress.  Among  these  are  closed  form  techniques  such  as 
a  Fredholm  or  convolution  integral,  and  approximate  methods  such  as  Runpe- 
Kutta.  Let  us  first  consider  a  first  order  Runge-Kutta,  sometimes  called 
Euler's  method.  In  this  method  the  differential  equation  is  solved  for  the 
stress  rate,  and  the  stress  terms  on  the  right  hand  side  are  assumed  at  time  t. 
Thus,  solving  equation  (29),  one  obtains 


da 

dt 


S 


t+dt  t  „  t+dt 

t+dt 


[If)  .  ^  + 


_  t+dt  t 

S  8° 


„  t+dt  t 

dSt  "  dSi  St  80 


dt 


dt 


/_3a_\  ga*" 

Wv . 


,  „  t+dt  t 

+  S  ga 


(dc  _  d£^  _  dcT\ 

ydt  dt  dt  J 


r„  t+dt 

S  8° 


■(*),  m . 


(jl\  gac  + 
L\3cP/ t 


_  t+dt  t 

St 


(36) 


I 
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where,  in  this  case 


t 

a 

n 


f  (at,et,Tt)  , 


(37) 


and  the  dashpot  coefficient  n  is  found  from  equation  (26)  and  (30)  to  be 


n  =  aC  g(ci,e,T)  . 


(38) 


Inspection  of  equation  (14)  reveals  that  if  one  defines 


1_  3F  3F 
t  9ct  3o  ’ 


(39) 


where  F  is  as  defined  in  equation  (11)  and  its  partial  derivative  is  evaluated 
at  a  =  a(t ),  then  equation  (39)  is  identical  to  the  plasticity  relation  given 
by  equation  (14).  To  complete  the  theory  it  is  necessary  only  to  assume  that 
g  and  f  are  zero  for  unloading.  The  impetus  for  using  a  third  order  expansion 
of  free  energy  in  conjunction  with  equations  (15),  (16),  and  (26),  as  well  as 
the  first  order  Runge-Kutta  approximation  is  therefore  justified  in  that 
introduction  of  these  assumptions  into  the  thermodynamic  constraints  [equa¬ 
tions  (4)  and  (5)]  will  bring  about  the  recovery  of  the  nonisothermal  plasticity 
theory  proposed  by  this  author.  The  conclusion  from  this  derivation  is  two¬ 
fold.  First,  the  classical  plasticity  theory  is  thermodynamically  consistent 
under  nonisothermal  conditions,  and  the  entropy  generation  term  a^  may  now  be 
defined  precisely  for  the  nonisothermal  plasticity  theory  proposed  by  the 
author.  Second,  the  nonlinear  viscoelasticity  theory  may  be  used  to  extend 
the  current  theory  to  model  more  complex  phenomena. 
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The  latter  statement  is  illustrated  by  the  fact  that  the  definition  of 


the  dashpot  coefficient  n  given  by  equation  (38)  requires  that 


Therefore,  the  recombination  of  the  plastic  and  creep  strain  rates  into 
a  single  inelastic  strain  produces  a  unified  theory  which  is  accommodated 
quite  simply  using  reduced  time. 

Now  consider  a  comparison  of  an  exact  solution  to  equation  (29),  as 
opposed  to  the  first  order  Runge-Kutta  approximation  given  by  the  plasticity 
theory  [equation  (14)].  For  this  purpose  we  introduce  reduced  time  [equation 
(24)]  into  the  differential  equation  to  obtain 


a  ,  1  da  _  de 

n  t+dt  dij;  “  d’l' 

*11 


t+dt  dij; 
fcM 


(41) 


Equation  (41)  may  be  Laplace  transformed  in  reduced  time  to  give 


i//  t+dt 


„  t+dt  r  >i  ,,  , ,,  ,,  t, 

a  =  EM  J  —  Of)  -  i]>  )  d (c  -  e  ) 
J  Hi!)* 


t+dt 

- 


+  cEt  f  <*-*')  dE« 

J n  xrr  d* 


where  if)'  is  the  dummy  variable  of  integration,  and  is  the  value  of  reduced 
time  at  time  t+dt.  Valanis  has  shown  that  a  single  Maxwell  element  is  in¬ 
sufficient  for  modelling  the  response  of  metals  [15].  However,  for  moderate 
load  increments  equation  (41)  is  adequate  since  the  generalized  Maxwell  model 
necessary  to  characterize  metal  response  consists  of  nonoverlapping  exponentials. 
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Thus,  a  single  exponential  is  predominant  at  any  strain  level,  with  the  appro¬ 
priate  Maxwell  element  described  by  the  current  material  data  in  equation  (42). 

The  strain  rate  and  elastic  modulus  rate  during  a  load  step  may  be  con¬ 
sidered  to  be  constant  inputs: 


d (c  -  rT) 
dt’ 


c  =  constant 


(43) 


and 


dt’ 


constant 


(44) 


so  that  equation  (42)  may  be  reduced  to  the  following  form 


t+dt 


t+dt 


t+dt •  f 

<”em  f  J 


t+dt 


dt'  +  r 


\f 


_em 


t+dt 


(\p-d’  ’ ) 


dt' 


(45) 


Therefore,  the  stress  increment  during  a  finite  load  step  is 


da  =  a  (t+dt)  -  o(t)  =  ^EMt+dtc+rE  ~  a^ 


4’  t+dt 

“N  (4-  -  *') 


e  n 


dvj' *  -  o(t)  .  (46) 


In  general  equation  (46)  must  be  solved  numerically  because  the  entropy 
generation  term  is  a  function  of  time.  It  should  be  noted  that  the  func¬ 
tional  form  of  the  nonlinear  viscoelastic  law  proposed  above  is  similar  to 
that  proposed  in  the  classical  plasticity  theory  [equation  (14)].  In  addition, 
its  application  to  a  finite  element  method  should  be  straightforward  due  to 
the  proper  rate  form. 
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A  COMPUTATIONAL  COMPARISON 


For  purposes  of  comparing  the  computational  efficiency  of  the  two  models 
an  example  problem  has  been  chosen  which  is  traditionally  quite  difficult  to 
solve  in  finite  element  codes.  The  problem  involves  a  material  not  unlike 
aluminum  with  isothermal  zero  time  stress-strain  curves  as  shown  in  Figure  6. 
At  elevated  temperatures  the  creep  strain  rate  may  be  approximated  by  a  power 
law: 


drC  =  9  x  1012  fo\ 5 
dt  T ( °K)  l  E  ) 


(47) 


The  presence  of  the  extremely  high  creep  rates  given  by  equation  (47)  requires 
that  most  finite  element  codes  solve  this  problem  using  very  small  thermo¬ 
mechanical  load  steps,  thus  consuming  considerable  computer  time.  This 
circumstance  is  due  to  the  fact  that  since  the  stress  increment  is  not  known 
a  priori,  the  creep  strain  rate  must  be  estimated  using  the  stress  at  the  start 
of  the  step. 

The  material  described  above  has  been  subjected  to  two  uniaxial  thermo¬ 
mechanical  load  histories,  as  shown  in  Figure  7.  Since  the  total  strain  for 
this  problem  can  be  determined  by  integrating  equation  (47)  and  adding  to  the 
elastic  and  plastic  strains,  the  solution  may  be  obtained  without  resort  to 
the  constitutive  laws  presented  herein.  However,  a  measure  of  the  accuracy 
of  the  constitutive  theory  is  to  apply  the  exact  total  strain  increment  to 
the  constitutive  theories  and  determine  which  better  predicts  the  proper 
stress  response.  This  problem  becomes  particularly  interesting  when  the 
thermomechanical  load  is  applied  at  slower  and  slower  rates,  thus  inducing 
rapidly  increasing  creep  strain. 


Under  the  conditions  noted  in  this  example,  if  one  assumes  that  the 
entropy  generation  term  is  constant  during  the  load  step,  which  is  identical 
to  plasticity  theory,  then  equation  (46)  may  be  integrated  in  closed  form  to 
give 


Valanis  [14]  has  succeeded  in  characterizing  the  reduced  time  <]'  for  isothermal 
loadings.  However,  since  the  plastic  strain  rate  [equation  (30)]  in  this  case 
is  both  temperature  and  stress  rate  dependent,  this  author  has  not  yet  found 
a  way  of  determining  this  parameter.  Therefore,  in  order  to  perform  this 
computational  comparison  it  was  necessary  to  assume  the  plastic  strain  rate 
in  equation  (26)  could  be  determined,  and  the  exact  values  were  utilized  in 
solving  equation  (48)  for  the  examples  presented  herein. 

In  addition,  two  other  methods  of  solution  were  tested  in  this  study. 
First,  a  fourth  order  Runge-Kutta  (Kut ta-Simpson  1/3  rule  [20]),  was  tested 
in  conjunction  with  equation  (14).  Second,  the  sub incrementation  approach  [21] 
was  used  with  equation  (14).  In  the  latter  method  one  cycles  through  the 
constitutive  law  several  times  using  sub incremental  loadings  on  each  load  step. 
For  these  tests  the  sub  increment  was  l  imited  to  a  total  strain  sub increment 
of  .0005  IN/IN  and  was  not  to  exceed  one  hundred  subincrements  during  the  load 
steps . 

Results  of  the  various  solution  procedures  are  shown  in  Figures  8  and  9, 
wherein  all  analyses  except  subincrementation  were  performed  in  a  single 
time  step.  The  results  were  obtained  by  first  applying  the  thermomechanical 
loading  in  the  time  period  shown  on  the  abscissa  and  determining  the  exart 
strain  increment  by  integrating  equation  (47)  and  adding  to  the  elastic  and 
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Figure  9.  Computational  Results  for  Test  Case  II 


plastic  increments  shown  in  Figure  6.  This  exact  strain  increment  was  then 
used  to  determine  the  predicted  stress  increment  in  each  of  the  four  models 
discussed  above.  All  of  the  models  are  seen  to  diverge  from  the  exact 
solution  for  decreasing  load  rate.  It  is  seen  that  the  fourth  order  Runge- 
Kutta  is  extremely  unstable  and  therefore  unusable.  This  is  due  to  the  fact 
that  the  differential  equation  becomes  increasingly  "stiff"  [22,  23]  for  in¬ 
creasing  time.  It  is  apparent  that  the  first  order  Runge-Kutta  and  nonlinear 
viscoelasticity  are  at  least  stable,  and  will  lead  to  an  over  approximation 
or  safe  estimate  of  the  stress.  The  subincrementation  method  is  demonstrated 
to  be  the  best  computational  scheme  for  the  examples  considered  here. 
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CONCLUSION 


The  author  has  compared  two  constitutive  theories  in  this  paper  on  the 
basis  of  physical  and  computational  accuracy.  It  is  seen  that  both  models 
satisfy  thermodynamic  as  well  as  physical  constraints.  In  addition,  both  may 
be  constructed  in  a  rate  form  suitable  for  application  to  finite  element  com¬ 
puter  codes.  Preliminary  evidence  indicates  that  subincrementation  used  with 
the  classical  plasticity  theory  better  satisfies  computational  requirements 
than  does  the  nonlinear  viscoelasticity.  However,  this  author  has  considered 
only  a  specific  type  of  nonlinear  viscoelastic  model,  and  final  -judgment  on 
the  issue  is  therefore  not  possible.  In  addition,  there  are  certain  aspects 
of  the  viscoelastic  model  such  as  the  single  integral  form  and  the  characteri¬ 
zation  of  state  variables  via  entropy  generation  terms  which  warrant  further 
studv. 
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